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Abstract. We derive a sharp, localized version of elliptic type gradient estimates for 
positive solutions (bounded or not) to the heat equation. These estimates are akin to 
the Cheng- Yau estimate for the Laplace equation and Hamilton's estimate for bounded 
solutions to the heat equation on compact manifolds. As applications, we generalize Yau's 
celebrated Liouville theorem for positive harmonic functions to positive eternal solutions 
of the heat equation, under certain growth condition. Surprisingly, this Liouville theorem 
for the heat equation does not hold even in R" without such a condition. We also prove a 
sharpened long time gradient estimate for the log of heat kernel on noncompact manifolds. 
This has been an open problem in view of the well known estimates in the compact, short 
time case. 



1. Introduction and main results 

In recent years, the following three types of gradient estimates have emerged to play a 
fundamental role in the study of the Laplace and the heat equation on manifolds. 

Theorem (Cheng-Yau [CY]). Let M be a complete manifold with dimension n > 2, 
Ricci(M.) > — k, k > 0. Suppose u is any positive harmonic function in a geodesic ball 
B(xq, R) C M. There holds 

(1.1) m<± + Cn ^ 

u R 

in B(xq, R/2), where c n depends only on the dimension n. 

Theorem (Li- Yau [LY]). Let M be a complete manifold with dimension n > 2, Ricci(M) > 
—k, k > 0. Suppose u is any positive solution to the heat equation in B(xq, R) x [to — T, to]. 
Then 

|Vu| 2 u* c n c n 

(1.2) J — 2 1 i <-^ + ^ + c n k, 

u A u R z 1 

in B(xq, R/2) x [to — T/2,to]. Here c n depends only on the dimension n. 

Theorem (Hamilton [H]). Let M be a compact manifold without boundary and with 
Ricci(M.) > —k, k > and u be a smooth positive solution of the heat equation with 
u< M for all (x,t) G M x (0, oo). Then 

, . |Vn| 2 ,1 , s , M 

1.3 J-jJ-< - + 2A; In— . 

u z t u 
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Clearly estimate (1.2) reduces to (1.1) when u is independent of time. On the other 
hand, for a time dependent solution of the heat equation, it is well known that the elliptic 
type gradient estimate (1.1) can not hold in general. This can be seen from the simple 
example where u(x,t) = e - '^ / 4 */(47rf) n / 2 being the fundamental solution of the heat 
equation in R™. The parabolic Harnack inequality also exhibits the same phenomenon 
in that the temperature at a given point in space time is controlled from the above by 
the temperature at a later time. However, in the case of compact manifolds, Hamilton's 
estimate (1.3) shows that one can compare the temperature of two different points at the 
same time provided the temperature is bounded. 

Just like the Cheng- Yau and Li-Yau estimates, it would be highly desirable to have a 
noncompact or localized version of Hamilton's estimate. However, the example in Remark 
1.1 below shows that the suspected noncompact version of Hamilton's estimate is false 
even for R™ ! This situation contrasts sharply with the Cheng- Yau and Li-Yau inequality 
for which the local and noncompact versions are readily available. 

Nevertheless in this paper we discover that, for noncompact manifolds the elliptic 
Cheng- Yau estimate actually holds for the heat equation, after inserting a necessary log- 
arithmic correction term. This correction term is slightly bigger than that in Hamilton's 
theorem (in the power of the log term). However the estimate holds for noncompact man- 
ifolds and it also has a localized version as the Cheng- Yau estimate. This result seems 
surprising since it enables the comparison of temperature distribution instantaneously, 
without any lag in time, even for noncompact manifolds, regardless of the boundary be- 
havior (see Remark 2.1). In some cases, our estimate (see (1.5) below) even holds for 
any positive solutions, bounded or not. This results seems new even in R™ or compact 
manifolds. 

Here is the statement of the theorem. 

Theorem 1.1. Let M be a Riemannian manifold with dimension n>2, Ricci(M.) > —k, 
k > 0. Suppose u is any positive solution to the heat equation in Qr,t = B(xq, R) x [to — 
T, to] C M x (—00,00). Suppose also u < M in Qr^t- Then there exists a dimensional 
constant c such that 

/-, ^ \Vu(x,t)\ ^ , 1 1 f-.. , M v 

(1.4) l - — / V < c(— + — rsr + Vfe)(l + ln— -) 

v ; u(x,t) - y R T 1 / 2 A u(x,t) J 

in Qr/2,t/2- 

Moreover, if M has nonnegative Ricci curvature and u is any positive solution of the 
heat equation on M x (0, 00), then there exist dimensional constants c±,C2 such that 

\Vu(x,t)\ 1 , , u(x,2t)v 

v ' u(x,t) - V/ 2V u(x,t) ' 

for all x G M and t > 0. 

An immediate application of the theorem is the following time-dependent Liouville the- 
orem, generalizing Yau's celebrated Liouville theorem for positive harmonic functions, 
which states that any positive harmonic function on a noncompact manifold with non- 
negative Ricci curvature is a constant. One tends to expect that Yau's Liouville theorem 
would still hold for positive eternal solutions to the heat equation. However the following 
simple example shows that this expectation is false. Let u = e x+t for x G R 1 . Clearly, u 
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is a positive eternal solution for the heat equation in R 1 and it is not a constant. Nev- 
ertheless, our next theorem shows that under certain growth conditions, Yau's Liouville 
theorem continues to hold for positive eternal solutions of the heat equation. Moreover, 
our growth condition in the spatial direction is sharp by the above example. 

Theorem 1.2. Let M be a complete, noncompact manifold with nonnegative Ricci cur- 
vature. Then the following conclusions hold. 

(a) . Let u be a positive eternal solution to the heat equation (i.e. solutions defined in all 

space-time) such that u(x,t) = e^^+V^*!) near infinity. Then u is a constant. 

(b) . Let u be an eternal solution to the heat equation such that u(x,t) = o(Jd(x) + y/\t\]) 
near infinity. Then u is a constant. 

Note that the growth condition on the second statement of Theorem 1.2 is also sharp 
in the spatial direction, due to the example u = x. 

Remark 1.1. Here we give an example showing that Theorem 1.1 is sharp for noncompact 
manifolds. This is suprising since it shows that Hamilton's estimate for the compact case 
(1.3) is actually false for noncompact manifolds. 

For a > consider u = e ax+a 1 . Clearly u is a positive solution of the heat equation in 
Q = [1,3] x [1,2] C R x (-00,00). Also Vit(2,2)/u(2,2) = a and M = sup Q u = e 3a+2a2 . 
Hence log(M/u(2, 2)) = a. Therefore at (x,t) = (2,2) the left hand side and right hand 
side of (1.4) with R = T = 1 are a and c(l + a) respectively. Obviously they are equivalent 
when a is large. 

Another application of Theorem 1.1 is the following sharpened estimate on the gradient 
of the fundamental solution of the heat equation. 

Theorem 1.3. Let M be a n dimensional complete, noncompact manifold with nonnega- 
tive Ricci curvature. Let G(x, y, t) be the fundamental solution of the heat equation. Then 
there exists c = c(n) > such that 

(lf?\ \V x G(x,y,t)\ 1 , | d(x,y) 2 ^ 

1 j G(x,y,t) - c t V* {+ t ] - 

holds for all x, y G M and t > 0. 

Remark 1.2. For a noncompact manifold and large time, it seems that all previous 
results establish an upper bound for |VG| rather than |VG|/G. Notice also that estimate 
(1.6) is almost sharp since, when M is the Euclidean space, one has 

\V x G(x,y,t)\ 1 d(x,y) 
G(x,y,t) 2yjt y/t 

We should mention that in the compact and short time case, estimate similar to (1.5) 
with a sharp power on the d(x, y) 2 jt term have been obtained for all derivatives of log G. 
See the papers [S], [H] (Corollary 1.3), [N], [MS], [Hs] and [ST]. The noncompact case has 
been wide open. 

Our work is motivated by [CY], [LY], [Yl-2], [H] and [BL] and [GGK]. 
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2. Proof of Theorems 1.1, 1.2 and 1.3 

The order of proofs is Theorem 1.1, 1.3 and 1.2. 



Proof of Theorem 1.1. 

Suppose u is a solution to the heat equation in the statement of the theorem in the 
parabolic cube Qr,t = B(xq,R) x [i — T, to]. It is clear that the gradient estimate in 
Theorem 1.1 is invariant under the scaling u — ► u/M. Therefore, we can and do assume 
that < u < 1. 

Write 

(2.1) / = ln«, ™=|Vln(l-/)| 2 = (|^2- 

Since u is a solution to the heat equation, simple calculation shows that 

(2.2) A/ + |V/| 2 -/ i = 0, 
We will derive an equation for w. First notice that 

_ 2V/(V/) t 2\Vf\ 2 f t 



(l-/)3 

2V/V(A/ + | V/| 2 ) + 2| V/| 2 (A/ + | V/| 2 ) 



(I-/) 2 (I-/) 3 
In local orthonormal system, this can be written as 



(2.3) Wt = JJJ i + 2 



Vj fiij + tfifjfij , J?fn + Wf\ A 
(I-/) 2 (I-/) 3 ' 



Here and below, we have adopted the convention uf = |Vu| 2 and u„ = An. 
Next 



( 2 - 4 ) Vw = (fl^Tyi) 

It follows that 



(l-/) 2 ^ (I-/) 2 (1-/) 



3 ' 



C9 ^ 2fifijj , ^fifijfj 

1 ' ' "(l-/) 2 + (l-/) 2 + (l-/)3 

^fifijfj / 4 2 /ii . R /? /j 



+ ,: ,v a +2 ,: tJ " +6 



(I-/) 3 (I-/) 3 (I-/) 4 ' 
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By (2.5) and (2.3), 

Aw — wt 

_ ^fij _j_ r jifijj ~ fjfiij 



(I-/) 2 (I"/) 2 

_l_ g IW| 4 _|_ g fifijfj _|_ 2 fifjj 



(1-fY (l-/)3 (1-/) 

fifijfj /j /jj r> 1^/1 



3 



4 _ 2- 



(l-/)3 (l-/)3- 

The 5th and 7th terms on the righthand side of this identity cancel each other. Also, by 
Bochner's identity 

fifijj ~ fjfiij = fjifjii ~ fiij) = Rijfifj — ~^\^f \ j 

where Rij is the Ricci curvature. Therefore 
Aw — wt 

(2.6) 2/g | g | V /|4 | g /■/■■/■ ^ /./../• o |V/| 4 2fc|V/| 2 



(I-/) 2 (I-/) 4 (I-/) 3 (I-/) 2 (I-/) 3 (I-/) 2 ' 
Notice from (2.4) that 

V fVw = + 2 - ^ 



(I-/) 2 (I-/) 3 " 
Hence 

(2.7) = 4 fifij {( - 2VfVw + 4- 



(I-/) 2 ' (I-/) 3 ' 
(2.8) = -4 /^^ + [2V/V W - 4- |V/|4 1 1 



(I-/) 3 L J (1-/) 3J 1-/- 
Adding (2.6) with (2.7) and (2.8), we deduce 
Aw — wt 

> + 2 |V/|4 + 4- /i/ ^ 



(l-/)3 

2 iVfl 4 2fe|Vf| 2 

^fij , l^-^l 4 1 /[ fifijfj > n 
(I"/) 2 (I"/) 4 (I"/) 3 " ' 

. 2/ |V/| 4 2fc|V/| 2 

Aw — w t > ^VfVw + 2 ir - 7F - Ir - 7F . 

Since / < 0, it follows that 

A -> 2 / V,T7 /, l V /| 4 M|Vjf 

A„ - „« > ^V/Vu, + 2(1 - njrw - W - J? , 



Since 



we have 



6 PHILIPPE SOUPLET AND QI S. ZHANG 

i.e 

2 f 

(2.9) Aw-w t > ^yV/Vw + 2(1 - f)w 2 - 2kw. 

From here, we will use the well known cut-off function by Li-Yau [LY], to derive the desired 
bounds. We caution the reader that the calculation is not the same as in that of [LY] due 
to the difference of the first order term. 

Let tp = ip(x, t) be a smooth cut-off function supported in Qr,t, satisfying the following 
properties 

(1) . V = ^(d(x,x ),t) = iP(r,t); ${x,t) = 1 in Q r/2 , t /a, < ^ < 1. 

(2) . V is decreasing as a radial function in the spatial variables. 

(3) . M < \m < ^when o < a < 1. 

(4) - |$ < ?■ 

Then, from (2.9) and a straight forward calculation, one has 
A(tfjw) + b ■ V(ipw) - 2^j- ■ V{ij)w) - (ipw) t 

(2 ' 10) iwi 2 

> 2^(1 - /)w 2 + (b ■ Vip)w - + (Aip)w - ip t w - 2kwip, 

where we have written 

2/ 

Suppose the maximum of ipw is reached at (xi, t\). By [LY], we can assume, without loss 
of generality that x\ is not in the cut-locus of M. Then at this point, one has, A(tpw) < 0, 
(ipw)t > and V(ipw) = 0. Therefore 

(2.11) 2^(1-/W(si,ti) < -[ (b-Vip)w -2^^-w + (A^)w -ip t w + 2kwip](x 1 ,t 1 ). 

W 

We need to find an upper bound for each term of the righthand side of (2.11). 

|(6 . wH < ^L|V/MW>| < 2w i l 2 \f \ |W| 



2,3/4 /IVVI 



[</,(! - /)]3/4 



This implies 

(2-12) |(6 • WH < (1 - f)i>w 2 + c R 4(f_ f) 3 - 

For the second term on the righthand side of (2.11), we proceed as follows 

IW| 2 /1/2 |Wf 

(2 13) ^ 

1 , 2 JVVf N 2 1,0 1 
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Furthermore, by the properties of i/j and the assumption of on the Ricci curvature, one 
has 

-(A^)w = -{d 2 ^ + (n - 1)-^ + 9 r ^9 r ln 

< (|d r 2 V| + 2(n - 1)^ + v^VI)™ 
<^J3^ + ^ B _ 1) JjW + , 1/V /|^ 

Therefore 

(2.14) -(Aif>)w< ^w 2 + c^ + ck^. 



Now we estimate \i/jt\ w. 



1,1 ,1/2 1^. 

|%| w = ip ' w 



^ 1 / / 1/2 \2 , / IV 7 * I \2 

This shows 

(2.15) |^<I^;2 +C _L 

Finally, for the last term, we have 

(2.16) 2kimj) < -4>w 2 +ck 2 . 

8 

Substituting (2.12)-(2.16) to the righthand side of (2.11), we deduce, 

2(1 - f)^w 2 < (1 - /)^ 2 + c RH l_ /)3 + + — + — + — + C fc 2 . 

Recall that / < 0, therefore the above implies 

/ 4 1 c c 

^W 2 (xi,tl) < C fi4 _ ^ 4 + -^ 2 (xi,tl) + + ^ + Cfc 2 . 

Since rrzjja < 1) the above shows, for all (x,t) in Q_r,t, 

il) 2 (x,t)w 2 (x,t) < ip 2 (xi,ti)w 2 (xi,ti) 
< ip(x 1 ,t 1 )w 2 (xi,t 1 ) 

C C ,9 

Notice that ip(x,t) = 1 in Qr/2,t/a and w = |V/| 2 /(1 - f) 2 . We finally have 



l-f(x,t)-R 

We have completed the proof of (1.4) since / = ln(u/M) with M scaled to 1. 
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To prove (1.5), we apply (1.4) on the cube Q^/i t / 2 = B(x,\/i) x [t/2,t]. By Li-Yau's 
inequality [LY], we know that 

M= sup u<cu(x,2t). 

Now (1.5) follows from (1.4). □ 
Proof of Theorem 1.3. 

Let G(x, y, t) be the fundamental solution of the heat equation on M. By the Li-Yau 
estimate [LY], given any 5 > 0, there exist ci, C2 > such that, 

(2.16) ^-p- e - W/M* < G{ t) < £i e -^,,)V(4+5)t 

V ' \B(x,Vt)\ ~ V ' " 

for all x,y G M and i > 0. Fixing x,y and t, we apply Theorem 1.1 on the function 
ix(z, r) = G(z, y, r) on the cube Q = B(x, s/i) x [t/2, t]. By the upper bound in (2.16) and 
the volume doubling property of the manifold, we know that 

( \ s Cl C3 



|S(z,v^)| - \B(x,Vt)\ 
for (2, r) G Q. By Theorem 1.1 and the lower bound in (2.16), we deduce 

< i (-> | 1n x 

U(x,t) - V/ 21 , , C 2 e -d( X , y )V(4-5)t)- 

\B(x,Vt)\ 

That is 

G(x,y,i) " C4 ^ ( + 



□ 



Proof of Theorem 1.2. 

(a). The proof is immediate. By our assumption the function u+1 satisfies ln(u + 1) = 
o(d(x) + s/i) near infinity. Fixing (xo,io) m space time and using Theorem 1.1 for u + 1 
on the cube B(xq, R) x [to — R 2 , to], we have 

|Vu(s ,td)| <C [1 + om 



u(x ,t ) + l R 

Letting R — > oo, it follows that |Vu(xo,io)| = 0. Since (xo,to) is arbitrary, one sees 
that u = c. 

(b). Fix (xo,to) in space time and let Ar = supn |it|. Consider the function U = 
■u + 2A2R. Clearly A 2 r < U(x,t) < 3A 2 r when (x,t) G Q 2 r^/ 2 r- Applying Theorem 1.1 
(a) for ?7, we deduce 

|Vu(x ,*o)| < £ 



u(x ,*o) +2A 2 r i?' 

Since A 2 r = o(R) by assumption, the result follows after taking R to 00. □ 

Remark 2.1. In case M has nonnegative Ricci curvature, (1.4) reduces to 

\Vu(x,t)\ . 1 1 . , , Ms 
^ c ^ + ^)( 1 + ln ^77T) 
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in Qr/2,t/2- If d(x,y) < \/T, we can integrate this along a geodesic to deduce 

u(y,t) y u{x,t)> 

for some 6 £ (0, 1) and c > depending only on n. Note in the compact case, Hamilton's 
theorem implies the above with = 1, which gives an elliptic Harnack inequality for the 
heat equation for compact manifolds. However, as pointed in the introduction, the elliptic 
Harnack inequality is false for the heat equation on noncompact manifolds in general. This 
explains the difference between our noncompact gradient estimate (1.4) and Hamilton's 
compact estimate (1.3). 
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